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Abstract 

The J-matrix inverse scattering nucleon-nucleon potentials (JISP), describing both two-nucleon 
£> ' data and bound and resonant states of light nuclei to high accuracy, are tested on the total photoab- 



o 



o 



X 



sorption cross sections of 2 H, 3 H and 3 ' 4 He. The calculations in the three- and four-body systems 



are carried out via the Lorentz integral transform method and the hyperspherical harmonics (HH) 



technique. To this end the HH formalism has been adapted to accommodate non-local potentials. 
The cross sections calculated with the JISP are compared to those obtained with more traditional 
realistic interactions, which include two- and three-nucleon forces. While the results of the two 
kinds of potential models do not differ significantly at lower energies, beyond the resonance peak 



they show fairly large discrepancies, which increase with the nuclear mass. We argue that these 
discrepancies may be due to a probably incorrect long range behavior of the JISP, since the one 
pion exchange is not manifestly implemented there. 
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I. INTRODUCTION 



In the last few years the possibility to perform very accurate calculations in few-nucleon 
systems has stimulated a new attitude in studying the nature and form of the nuclear force. 
Differently from what was common in the past, nowadays potential models are tested not 
only on the increasingly accurate nucleon-nucleon (NN) data and triton binding energy, but 
also on bound and resonant states of more complex systems. 

The modern debate on nuclear forces focuses on potential models that can schematically 
be grouped into three different categories. Category A includes the most traditional poten- 
tials. Their two-body parts are either based on meson exchange models Q, E , [J , or largely 
phenomenological |5| . However, all of them include the essential one pion exchange term 
at long range and reproduce the NN scattering data with very high accuracy. However, 
when tested on A > 3 nuclei, they show non negligible underbinding. Three-body forces 
have been introduced to fit the A=3 binding energy, and more recently also to describe 
energies of bound and low-lying resonant states of p-shell nuclei jfj Q, si 0]- They require 
very intense computational efforts when used in may-body systems. 

Category B includes the effective field theory (EFT) potentials Q], which have a stronger 
connection to Quantum Chromo-Dynamics (QCD). They imply consistent two-, three- or 
more-body terms. The number of these terms and of the associated parameters, as well as 
the accuracy in reproducing data, increase significantly with the order of the expansions. 
The common feature of potentials belonging to categories A and B is the presence of three- 
body operators. These require very intense computational efforts when used in many-body 
systems. 

Category C includes potentials that, while describing NN scattering data with high accu- 
racy make use of the remaining off-shell freedom in order to reproduce few-body ground 
state e nerg ies and low-lying resonances, avoiding the necessity of three-body operators 

HQffUQ. 

Electromagnetic reactions can provide good constraints of the off-shell behavior of the 
potential. The potentials of category A have been tested at large on quite a number of 
observables in the photo- and electro-disintegration of the two-body system. The comparison 
with data is very good, at least within the non relativistic regime (see e.g. As 
to category B potentials, analogous systematic studies are still missing. Investigations on 
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electromagnetic reactions have been performed also in the three-body systems (see e.g. 



18, 
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20 



21 1) for potentials falling into category A and B. 



Due to its large binding energy the four-body system represents a very good testing 
ground for off-shell versus many-body properties of the nuclear force. However, only re- 
cently an ab initio calculation of the total photodisintegration cross section of 4 He has been 
performed 22 1 with a category A potential of local nature (AV18 and UIX) P, ^|. Even if 



the lack of precise experimental data limits the validity of such kind of calculations as test 
of the off-shell behavior of the potentials, it is interesting to calculate the same observable 
also with a potential of category C. In particular, one can ask whether potential models that 
reproduce the binding energy of the alpha particle in two different ways, either by includ- 
ing three-body operators, or by a proper choice of the off-shell behavior of the two-body 
potential, predict the same results for the electromagnetic observables. 

Among the potentials of category C, of particular interest are the potentials of Refs. [3, 

Q| , which are highly non-local and given in terms of matrix elements constructed by 
means of the J-matrix version of inverse scattering theory. Due to the lack of the three-body 
operator and the way they are constructed, they lead to rapid convergence, when used in 
many-body ab initio calculations based on complete basis expansions, as was shown in k| 
within the no core shell model approach 23^ . 

The aim of this work is twofold. On the one hand we want to extend the calculation of 
Ref. 22] to non-local potentials. On the other hand we would like to compare the results on 
the total photodisintegration of 2 H, 3 H and 3 ' 4 He below pion threshold, obtained with two 
versions of the JISP, to those given by typical category A potentials. 

The choice of this particular electromagnetic observable is that in this case the validity of 
Siegert's theorem allows to take largely into account, in an implicit way, the two- (or more-) 
body current contributions. Because of the particular form of the JISP (in matrix and not 
in operator form) it would be otherwise very difficult, if not impossible, to construct them 
in a consistent way. On the other hand, just these currents can help to get an idea on the 
importance of the underlying degrees of freedom that are "mocked" by the off-shell part of 
the potential. 

The photodisintegration cross sections are calculated here using the Lorentz Integral 
Transform (LIT) method 2m which allows the treatment of the continuum states dynamics 
by means of a bound state technique. This method is briefly summarized in Sec. II. We use 
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the symmetrized HH expansion 

|2a I2J, 27, 28] for the three- and four-body calculations, 
while standard numerical methods are used to solve the differential equations in the two- 
body case. 

The use of the HH formalism is natural in configuration space and for local potentials. 
Since the JISP are non-local and given in terms of two-body matrix elements between 
harmonic oscillator (HO) states, the many-body HH formalism has to be adapted to this 
case. Sec. Ill contains the description of how this is achieved. In Sec. IV we present 
the results for the total photodisintegration of A=2,3,4 nuclei and compare them to those 
obtained by category A potentials. Conclusions are drawn in Sec. V. 

II. THE TOTAL PHOTOABSORPTION CROSS SECTION WITH THE LIT 
METHOD 

The total photoabsorption cross section is given by 

a 7 (cj) = 4ir 2 aujR(uj) , (1) 

where a is the fine structure constant, u represents the energy transferred by the photon 
and R(u) is the response function defined as 

^) = il(/|0|O)| 2 ^/-^-^o). (2) 



Here |0) and Eq are the nuclear ground state wave function and energy, |/) and Ej denote 
eigenstates and eigenvalues of the nuclear Hamiltonian H, and is the operator relevant 
to this reaction (it is assumed that recoil effects are negligible). In the low-energy region 
considered here one can rely on Siegert's theorem and use for the unretarded dipole 
operator 



e^E^f- (3) 



A 

2 

The calculation of R(u) seems to require the knowledge of the continuum states \n). 
However, this difficulty can be avoided usin g t he LIT method. This method has been 
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described extensively in several publications 
steps which are needed for the calculation of R{ui). 
Step 1. The equation 



Here we only summarize the three 



(H-E o -u o + ir)\$) = e\0) (4) 



has to be solved for many ujq and a fixed T. This is a Schrodinger-like equation with a 
source. It can be shown easily that the solution is localized. Thus one only needs a 
bound state technique to calculate it. One generally adopts the same bound state technique 
as for the solution of the ground state, which is an input for Eq. We use expansions on 
the HH basis, as explained in Sec. III. 

The values of the parameters ujq and T are chosen in relation to the physical problem. 
In fact, as it becomes clear in Step 2, the value of T is a kind of "energy resolution" for 
the response function and the values of Uq scan the region of interest. In our case we are 
interested in the resonance region and up to pion threshold. Therefore we solve Eq. (0J) with 
r = 10 and 20 MeV and for a few hundred of Uq values chosen in the interval from —10 
MeV to 200 MeV. 

Step 2. After solving Eq. (j2J) the overlap (^l^) is calculated. Of course this overlap 
depends on uj and T. A theorem for integral transforms based on the closure property of 
the Hamiltonian eigenstates (^J ensures that this dependence can be expressed as 24] 

L(cj , r) = ) = J R{uj) C(u, u , T) duo , (5) 

where C is the Lorentzian function centered at ujq and with T as a width: 

C(u,u ,r) = 1 . (6) 

(u - u y + T 2 



Therefore by solving Eq. (j3J) one can easily obtain the LIT of t 



Step 3. The transform (jHJ) is inverted (see e.g. Ref. 0, 



re response function. 



8ll |) in order to obtain the 
response function and therefore the cross section. Of course the inversion result has to be 
independent of V and show a high degree of stability. 



III. ACCOMMODATION OF NON-LOCAL NN POTENTIAL MATRIX ELE- 
MENTS IN THE HH FORMALISM 

In this section we describe the implementation of non-local potentials in the HH formal- 
ism. Our method is alternative to that of Ref. |32j and more convenient for the structure of 
our codes. 
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A. The A-body basis states 



In the HH formalism the antisymmetric A-body configuration-spin-isospin basis func- 
tions with total angular momentum J A , J A and isospin T A ,T A are given by (notations as in 
Ref. 



\n A K A J A J z A T A T z A T A a A f3 A ) 
At 



n A >^^=^ \K A L A M A T A Y A ^a A )\S A S\T A T z A T A ,Y A ^p A ) " A , (7) 



Ja 



J J A 



where 

(f^A I K A L A M A Y A Y A ^ia A ) = y^M^v^a^A) (8) 

are the HH functions with hyperspherical angular momentum K A , and orbital angular mo- 
mentum quantum numbers L A , M A that belong to well defined irreducible representations 
(irreps) T 1 E T 2 ■ ■ ■ E T A of the permutation group-subgroup chain S% C «S 2 . . . C S A , de- 
noted by the Yamanouchi symbol [T A , Ya-i] = [Fa, Ta-i, . . . , Ti]. The dimension of the 
irrep T m is denoted by |T m | and A rAj y A1 is a phase factor 3^. Similarly, the functions 

(sI-.s^.^SaS^TIY a ,Y a _ x (3 a ) ee X^ wiFa ^_ i/3a (^..4^i--^) (9) 

are the symmetrized spin-isospin basis functions. In the definitions (JBJ) and Q the quan- 
tum numbers a A) (3 A are used to remove the degeneracy of the HH and spin-isospin states, 
respectively. The function 

(pM = R [ Ap) (10) 

is the A-body hyperradial basis function, 



where L^(x) are the associated Laguerre polynomials. The A-body hyperradial basis func- 
tions depend on the range parameter b and the Laguerre parameter ocl- 



For the sake of brevity in the following the state of Eq. will be denoted by \n A , [K A ]), 
with [K A ] = K A J A J\T A T z A T A a A (3 A . 
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B. The interaction and the HH (A — 2, 2) basis 



The representation of a non-local two-body interaction becomes very simple if the basis 
for the v4-body Hilbert space is chosen as an outer product of a 2-particle and a (A — 2)- 
particle states. In such a representation the interaction has the following form 



v [ %3)= £ \c2Ca-2)v [ ^ 2 (c 2 C a ^\, (12) 

C 2 c' 2 Ca-2 



where by c 2 and c' 2 we denote the two-body states of the interacting particles ij and by 
Ca-2 we indicate a state that includes the (A — 2)-body residual system. For a fermion 
system, these states should be taken as antisymmetric two-body and (A — 2)-body states, 
respectively. 

We choose to label the particles of the pair as particle A and A — 1 and to indicate their 
relative coordinate by ff = f)A-i = J\{ta—Ta-x)- Thus c 2 depends on the Jacobi coordinate 
774-1 and Ca-2 on the remaining (A — 2) reversed order Jacobi coordinates tJa-2,Va-3---Vi 
|27| . Notice that in this way Ca-2 has information also on the relative orientation of the 2- 
and the (A — 2)-body subsystems. 

After having transformed {?7a-2> Va-3 ■■■ Vi} and V to hyperspherical and spherical coor- 
dinates, respectively, one can express Ca-2 (C2) on the hyperspherical (spherical) basis. In 
analogy to Eq. (J2J) one has the antisymmetric (A — 2)-body configuration-spin-isospin states 

Ar 



Ca-2) = \riA-2) ^2 



^F A . 2 ,Ya-3 



y A _ 3 V\ t a- 2 \ 

~ ~ 1 Ja-2 

\K A -2L A-2M a-2 1 ? A-2Y A-Z& A-2) \S A-^A-iF A-2Xa-2^ A-2-, Y A -Z&A-2) z , (13) 

and the antisymmetric two-body states (the sum £ 2 + 5*2 + T 2 is odd) 

|c 2 ) = \n 2 )[\t 2 m 2 )\S 2 S z 2 T 2 T>)]% . (14) 

The notations in the previous equations are similar to those in Eqs. (J7J) and (8). The 
difference is in the subscripts A — 2 and 2, indicating that the quantum numbers refer to 
the (A — 2) system and to the pair, respectively. Therefore the corresponding HH, spin- 
isospin and hyperradial functions are ~y\t A _} 2 LA-2MA-2TA-2YA-3 a A 

(Qa-2) and Yi 2m2 (Q 2 = fj), 

x' A ~ 2] 9Z rp TZ fi ?a Ba {s\..s z A _ 2 ,t\..t\_ 2 ) and X [ s 2 ssT 2T A s i4,ttt%), i^l"? (pA-2) and 

b A-2 b A-2 1 A-2 1 A-2 1 A-2* A-3PA-2 ^2^2 

RjliJjl), respectively. Moreover, in the following \Ca-2) will be denoted by \ua-2, [Ka-2\) 
and |c 2 ) by |n 2 , [^2]}- 



In order to calculate the matrix elements of the potential between the A-body states 
of Eq. (J7|), it is also useful to introduce a particular HH basis that reflects the division 
of the particle system into a pair of interacting particles and an (A — 2) spectator. Such 
a basis is obtained by coupling the 2- and (A — 2)-body spherical-spin-isospin states and 
hyperspherical-spin-isospin states, to yield an A-body basis function with quantum numbers 
[K A - 2 ], [K 2 ],Ka, Ja, J a, T a , T|, through the relation 

\[K A -2}; {K 2 ])K a JaJ z a TaTX) = AC >b (sin^ A _ 2)2] )^(cos^^ 2 , 2] )^- 2 Pi a - b )(cos2^ A _ 2)2] ) 

x[\[K A - 2 ])\[K 2 ])]%%, (15) 



where 6[a-2,2] is defined by the following equations 



7] = p Sin 0[A- 2 ,2] 
PA-2 = pCOs9 [A -2,2} ■ (16) 



Here Pi a ' b) are the Jacobi polynomials, with arguments 

a = £ 2 + 1/2 

b = K A -2 + — - — 

K A - K A -2 - ?2 

n = . 

2 

The numerical factor 



(17) 



b _ 2(2n + a + b)n\T(n + a + b+l) 
n V T(n + a + l)T(n + b + l) U8J 

is a normalization constant. 

For the sake of brevity we use the following notation: 

P^- 2/2 = AC' fe (sin^ [A _ 2i2] ) fe (cos0 [A _ 2 , 2] )^- 2 P^(cos20 [A _ 2 , 2] ) . 
Therefore the HH (A — 2, 2) basis states |[i^ A _ 2 ]; \K 2 \)K aJ aJ z a^ a^aI can now be written as 



\([K A -2}; [K 2 })K a JaJ a TaT z a ) = P^ Wa-2]) \[K2])] J jf T i • (19) 

These states, together with the A-body hyperradial basis states |n A ), form a complete or- 
thonormal basis of our Hilbert space. 



S 



C. The transformation between the A and (A — 2,2) basis 



The next step is the evaluation of the overlaps between the A-body functions, Eq. ((7J), 
and the (A — 2)— and two-body functions, Eq. (J13j) and (|14|) . To this end we use the 
completeness of the HH (A — 2, 2) basis, Eq. (fTU|) . i.e. 

(n A , [K A \ I n A - 2 , [K A _ 2 ]; n 2 , [K 2 \) = ^(tia, [K a ] \n A , ([K A - 2 ]; \K 2 })K A J A J Z A T A T Z A ) 

x (n A ([K A - 2 ]] [K 2 ])K A J A J Z A T A T Z A | n A _ 2 , [K A _ 2 ]; n 2 , [K 2 \) . (20) 

Let us start with the first matrix element on the right hand side of Eq. (J2(Jj) . The contribu- 
tions of the hyperangular, spin and isospin matrix elements are evaluated with the help of 
6j and 9j symbols, the hyperspherical coefficients of fractional parentage (CFPs) and the 
spm-isospm CFPs 133. One has 



([K A }\ {[K A - 2 ]) [K 2 ]) K A J A J A T A T A ) 



^{2J A _ 2 + 1)(2J 2 + 1){2S A + 1)(2L A + 1) < 



s 2 S A - 2 S A 
h L A _ 2 L A 
■h J A -2 J A 



^ ^TaFa-Ava-iFa- 



x (K A L A Y A a A \ (K A _ 2 L A _ 2 Y A _ 2 a A ^ 2 ; £ 2 ) K A L A ) 
x {S A T A Y A (3 A \ S A _ 2 T A _ 2 Y A _ 2 (3 A _ 2 ; S 2 T 2 ) S A T A ) 



(21) 

where the hyperspherical matrix elements are written as 

{K A L A Y A a A \ {K A _ 2 L A _ 2 Y A _ 2 a A ^ 2 ; l 2 ) K A L A ) = 

(K A L A T A _ia A ^i\K A L A T A a A ) 

x ((K A _ 2 L A _ 2 T A _ 2 a A - 2 ; t 2 )K A L A \K A L A Y A „ x (x A _ x ) (22) 

and, analogously, the spin-isospin term is written as 



(S A T A Y A f3 A | [S A ^ 2 T A _ 2 Y A _ 2 (3 A ^ 2 ; S 2 T 2 j S A T A ) — 

^ n tscfp a 

Sa-iTa-iPa-i 

x {S A - 2 S A -iS A \{S A - 2 \ S 2 )S A )(T A - 2 T A -iT A \(T A _ 2 ; T 2 )T A ) 

with CFPs products 

T1TSCFP A = (SaSa-iTaTa-xTa.xPa^ISaTaTaPa) 



(23) 



x (Sa~iSa-2Ta~iTa-2^a-2Pa-2\Sa-iTa-i^a-iPa-i) ■ (24) 

The spin-isospin matrix elements can be easily evaluated using angular momentum tech- 
niques: 

(T A -2Ta-iTa\(Ta-2;T 2 )Ta) = (-) Ta+t --> + W(2Ta-i + 1)(2T 2 + 1) \ ^ * Ta_1 1 . 

[ * T A T 2 j 

(25) 

Here i is equal to | and stands for the isospin of a single nucleon. The spin matrix element can 
be obtained by simply replacing in Eq. (|25|) the isospin quantum numbers by corresponding 
spin quantum numbers. The overlap between the interaction basis and the HH (A — 2, 2) 
states are given by 

(n A {[K A -2]\ [K 2 ])K A J A J Z A T A T Z A I n A -2, [K A - 2 ]; n 2 , [K 2 ]) = 

(Ja-2J%^J2J z 2 \JaJ%)(T A -2TX^T 2 T^\T a ^ (26) 

where, 

A n n ^ n *(K A ,K A _ 2 J 2 ) = 

J v 2 dvR®e 2 (v) J P A A S 2 7 d P A-2 4t?(^ 2 )^1(P) P ?r 2A ( c ^[^ 2 ,2]) • (27) 

This two dimensional integral can be easily evaluated using the Gauss-quadrature integra- 
tion. 



D. The potential matrix elements 

Using the representation of the potential as in Eq. (jl2j) and the overlaps 
(n A , [K A ] | c 2 C A - 2 ) calculated above, the A-body matrix elements of a scalar-isoscalar 
two-body operator can be written as 

(n A [K A ]\v®(A, A - 1) | n' A [K' A \) = E (n A [K A ) \c 2 C A „ 2 )v^ c , (c' 2 C A . 2 \ n' A [K' A ]) = 

C 2 c' 2 Ca-2 

E E KV n2 ( K A,K A . 2 ^ 2 )^ n \K' A ,K A ^e 2 ) 

[K 2 ],[K 2 ],[K A _ 2 ] [n 2 ],[n' 2 ],[n A _ 2 ] 

x ([K A ] | {K A _ 2 \- [K 2 })K A J A J Z A T A T Z A )([K' A ] \ ([i^_ 2 ]; \K' 2 ])K' A J A J Z A T A T Z A ) 

<[* 2] ,W^-1), (28) 
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where we have used the orthogonality of the Clebsch Gordon coefficients 

/] {Ja-2Ja-2^ J 2\JaJa){Ja-2Ja-2 J 2J2\J'a J a) = $J A J'Jj%J'I ( 29 ) 

JZ JZ 

J 2 J A-2 

and 

E/rri rjlZ rji rriz I rri rriz \ /rri rjlZ rri rfiz I rjl/ rjilz \ c £ / OC\\ 

{lA-2±A-2 1 2J-2 \ 1 A±a)\ 1 A-2J-A-2 1 2J-2\ 1 A 1 a) = °T A T^O T z T u . (dU) 

1 2 

Finally, the actual potential matrix elements are calculated via 



v 



[2] 

n 2 [K 2 },n' 2 [K' 2 ] 



h 2 r^j^T 2 vMrAn^ (31) 



E. Application to the J-matrix inverse scattering potentials 

In view of the use of the JISP as given in 

HQ 

, a convenient choice for the configuration 

space two-body basis functions are the HO states, 

(r> 2 £ 2 m 2 >HO = Rn% 2 {r)Y l2m2 {r) , (32) 

with the radial function 

- ^f^Mrm Q'W-v^ieV/rS) . (33) 

Here r is the distance between the two particles i.e. r = \ fi— f^l, ro is related to the oscillator 
strength VL HO an d the reduced mass /i through the relation r = \Jh/ [iVLho- However, the 
basis states in Eq. (|3*2*|) do not possess the proper normalization. In fact they are normalized 
according to 

J df\(f]nlm) HO \ 2 = l, (34) 

while in the canonical transformation from the single particle coordinates to the centre of 
mass and reversed order Jacobi coordinates the appropriate normalization is 



drf\ (ff\nlm) | =1, (35) 

where fj = \I\{ta — Ta-i) = \J\t- Consequently a normalization factor of v^8 has to be 
introduced. Therefore the appropriate two-body basis functions are 

(V\n 2 £ 2 m 2 ) = R®MYe 2m M ee v^<£(r = v 7 ^) W?) • (36) 
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IV. RESULTS 



In the following we discuss the results for the unretarded total photoabsorption cross 
section of the two-, three-, and four-nucleon systems obtained with the JISP and compare 
them to the results given by other nuclear force models. For the three- and four-body 
systems the cross sections are calculated with the LIT method, as outlined in Sec. II, and 
using HH expansions of ground state wave function and For the JISP these HH expansions 
have rather rapid convergences, though no effective interaction has been introduced in this 
case. We obtain convergent binding energies of the three- and four-body systems with 
K < 20 (see table EJ), while for potentials like AV18 one needs K larger than 100. For the 
deuteron photoabsorption with the JISP we also use the LIT method, but here combined 
with expansions on HO functions. For the other NN potentials a conventional calculation 
with explicit NN final state wave functions is carried out (for one of these potentials it was 
checked that we get the same result also with the LIT method). 

Before discussing the total photoabsorption cross sections we first consider the longitudi- 



121, AV14 



nal deuteron elastic form factor Ft (a < ). In Fig. 1 we show results for the JISP6 

n n ~ 

|4J], and BonnRA [1| potential models (note that the proton and neutron electric form factors 
are set equal to 1 and 0, respectively). Up to a momentum transfer q 2 = 2 fm -2 all models 
lead to rather similar results, while for higher q 2 the JISP6 result is considerably lower than 
those of the other two potential models. In comparison with experimental data one usually 
considers the elastic deuteron form factors A(q 2 ) and B(q 2 ), where A(q 2 ) contains besides 
the longitudinal also a transverse contribution. However in the momentum range shown 
in Fig. 1 the transverse piece is very small (see, e.g. For this momentum range it 

is known that the A(q 2 ) of conventional potential models like AV14 or BonnRA agree well 
with experimental data (see e.g. jj]). Thus one may conclude that the JISP6 potential leads 
to a good description of data only up to q 2 = 2 fm -2 . On the other hand this might be suf- 
ficient for the calculation of electromagnetic low-momentum transfer reactions like nuclear 
photoabsorption below pion threshold. 

In Fig. 2 the unretarded deuteron total photoabsorption cross section is shown for various 
potential models. In Fig. 2a one sees that the low-energy peaks have slightly different 
heights (note that the cross section below about 5 MeV is also affected by the here not 
considered Ml contribution). Beyond the peak one finds a good agreement of all models 
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up to about 10 MeV. At higher photon energies one has very similar results for the AV14, 
rn 

AV18 |5j, BonnRA models, while the JISP6 result shows a somewhat higher cross section 
between 15 and 60 MeV (see Fig. 2b). Fig. 2c helps to understand the situation better. 
For the BonnRA potential we illustrate the effect of the inclusion of retardation and of 
other multipoles (electric and magnetic ones up to multipole order L—4), where the non- 
relativistic one-body current and the Siegert operator are considered (calculation corresponds 



to Arenhovel's normal calculation for the deuteron photodisintegration, see e.g. [37[). One 
notes that the additional contributions increase the cross section only slightly. In principle 
one has to include also other current contributions (meson exchange, isobar, and relativistic 
currents), however, as shown in Fig. 7.1.8 of Ref. |37j these effects are rather small. In the 
considered energy range they lead to a reduction of the normal cross section between 2% and 
3%, thus bringing the total result very close to the unretarded total photoabsorption cross 



section. In comparison to experimental data |38L I39L I40L l4ll l42j ] there is good agreement of 



the unretarded cross section with the BonnRA potential, while with the JISP6 potential the 
data are somewhat overestimated . 

Now we turn to the total photoabsorption cross section of 3 H and 4 He. In Fig. 3a we 
show the triton case with the JISP6 and JISP16 potentials and for the nuclear force models 



and Tucson-Melbourne 



M) 
4c 



AV14 Q and UVIII Q; AV18 ^ and UIX J BonnRA Q 

0. Results for the AV14+UVIII and BonnRA+TM potentials are taken from Ref. 
The situation is similar to the deuteron case of Fig. 2: small differences among the low- 
energy peak heights of all potential models and beyond 20 MeV differences between the 
conventional nuclear force models and the JISP, which, however, are considerably more 
pronounced than for the deuteron case. Since beyond the peak there are no triton total 
photoabsorption cross section data available, the comparison with data is done for the 3 He 
case (see Fig. 3b), choos,ng one of the JISP (JISP6) and one eonvent.ona. (AViS+UIX) 
potential model. Unfortunately the data 44] have rather large error bars and do not lead 
to a clear picture. However, one can say that between 20 and 50 MeV there is a better 
agreement with the JISP6, while beyond 50 MeV the conventional potential models are 
favored. 

In Fig. 4 we show the 4 He total photoabsorption cross section for the JISP6 and JISP 16 



potential mode 
is described in 



s, as well as for the AV18+UIX nuclear force. The calculation for AV18+UIX 
22l | . There it is discussed that the convergence of the HH expansion is rather 
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slow because of the three-nucleon force. As already mentioned, for the JISP we find a much 
better convergence, though no effective interaction is introduced, in Fig. 4 one notices that 
the trend observed for deuteron and three-nucleon cases is confirmed. Beyond the peak the 
cross sections obtained with the JISP are considerably higher than that obtained using the 
conventional nuclear force model AV18+UIX. This effect is more pronounced than for the 
three-nucleon case, which in turn is already stronger than found for the deuteron case. In 
addition, for 4 He the JISP lead to somewhat lower peak heights and to a shift of the peak 
position by about 2 MeV towards higher energy, when compared to the AV18+UIX result. 



It is evident that the latter agrees much better with the experimental data 4j|. Here we 
should mention that the experimental situation is not yet settled (see also 221]). particularly 
regarding the height of the low-energy peak (the shown data are the only experimental total 
He a 1 (uj) results that extend to energies beyond the peak region). 



i 



V. CONCLUSIONS 



In this work we have presented results for the total photoabsorption cross section of the 
two- three- and four-nucleon systems, obtained within the recently proposed J-matrix inverse 
scattering potential models. 

The calculation in the three- and four-body systems are performed via the LIT method 
and the HH technique. To this aim the HH formalism has been adapted to accommodate 
non-local potentials. 

The comparison of the JISP results with those obtained with a few more traditional 
realistic potentials, including two- and three-body forces, has shown that, while the latter 
give very similar results, JISP display a rather different behavior, especially beyond the 
resonance peak. The differences increase with the nuclear mass. 

These results show that a " classical" electromagnetic observable as the total photonuclear 
cross section is able to emphasize the rather different off-shell properties of these two classes 
of potentials. In particular, in this observable non-localities of the two-body potential and 
three-body forces are not equivalent. Considering that, different from the JISP, the tradi- 
tional potentials have in common the long range one-pion exchange part, one may ascribe 
the discrepancies to this fact. 

In the case of deuteron the comparison with data speaks in favor of the more traditional 
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two-body potentials. For the three-body systems data are missing or have insufficient accu- 
racy, making the comparison inconclusive. Regarding the four-body system the old data of 
Ref. H seem to favor the traditional potentials, as in the deuteron case. More recent data 
obtained in Lund 4||, while measuring only the 4 He(7, n) cross section point in this direc- 
tion as well, provided that they are extrapolated to <r 7 by means of the charge symmetry 
argument. 

More data of higher accuracy are very much needed to further clarify the issue of two- 
and three-body forces. 
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FIG. 1: Deuteron elastic longitudinal form factor (see text) for three different potential models. 
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Fig. 2 




FIG. 2: Deuteron total photoabsorption cross section at lower (a) and higher (b) photon energies 



for four different potential models. Comparison with data for two potentials (see text) is shown in 
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Fig. 3 
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FIG. 3: Three-body total photoabsorption cross sections: (a) 3 H results for five different potential 
models; (b): 3 He results for two potential models, compared to data from Ref. \^4\. 
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TABLE I: Binding energies [MeV] for A = 3, 4 nuclei using the HH expansion with the JISP6 and 
JISP16 potentials. For comparison we also present the results obtained in the NCSM approach. 



Nucleus 


JISP6 
HH NCSM [12] 


JISP16 
HH NCSM [14] 


3 H 


8.461(1) 8.461(5) 


8.369(1) 


8.354 


3 He 


7.749(1) 7.751(3) 


7.662(1) 


7.648 


4 He 


28.602(1) 28.611(41) 


28.299(1) 


28.297 
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